Mon. Not. R. Astron. Soc. 000, 1-21 (2011) Printed 11 November 2011 



(MN BTeX style file v2.2) 



Spatial matter density mapping of the STAGES Abell 
A901/2 supercluster field with 3D lensing 



4 



P. Simoni'2, C. Heymans^, T. Schrabback^, A.N. Taylor^, M.E. Gray 
L. van Waerbeke^ C. Wolf^, D. Bacon^ M. Barden^ A. Bohm^ B. Haufiler^ 
K. Jahnke^, S. Jogee^^, E. van Kampen^^, K. Meisenheimer^, C.Y. Peng^^ 

^ Argelander-Institut fiir Astronomie, Universitdt Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany. 

^ The Scottish Universities Physics Alliance (SUPA), Institute for Astronomy, University of Edinburgh, Edinburgh EH9 3HJ, UK. 

'"^Leiden Observatory, Leiden University, Niels Bohrweg 2, NL-2333 CA Leiden, The Netherlands. 

"^School of Physics and Astronomy, The University of Nottingham, University Park, Nottingham NG7 2RD, UK. 

^Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, V6T IZl, Canada. 

^Department of Physics, Denys Wilkinson Building, University of Oxford, Keble Road, Oxford, 0X1 3RH, UK. 

"^Institute of Cosmology and Gravitation, University of Portsmouth, Dennis Sciama Building, Burnaby Road, Portsmouth, POl 3FX, UK. 
^ Institute for Astro- and Particle Physics, University of Innsbruck, Technikerstr. 25/8, A-6020 Innsbruck, Austria. 
Max- Planck- Institut fiir Astronomic, Konigstuhl 17, D-69117, Heidelberg, Germany. 

University of Texas, McDonald Observatory, Fort Davis, TX 79734, USA. 

European Southern Observatory, Kari-Schwarzschild-Str 2, D-85748 Garching, Germany. 

NRC Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria, V9E 2E7, Canada. 



ACCEPTED SEPTEMBER 2011 



ABSTRACT 

We present weak lensing data from the HST/STAGES survey to study the three- 
dimensional spatial distribution of matter and galaxies in the Abell 901/902 superclus- 
ter complex. Our method improves over the existing 3D lensing mapping techniques 
by calibrating and removing redshift bias and accounting for the effects of the radial 
elongation of 3D structures. We also include the first detailed noise analysis of a 3D 
lensing map, showing that even with deep HST quality data, only the most massive 
structures, for example M200 ^ lO^^MQh"^ at z ~ 0.8, can be resolved in 3D with any 
reasonable redshift accuracy (Az « 0.15). We compare the lensing map to the stellar 
mass distribution and find luminous counterparts for all mass peaks detected with a 
peak significance > 3(t. We see structures in and behind the z = 0.165 foreground 
supercluster, finding structure directly behind the A901b cluster at z ~ 0.6 and also 
behind the SW group at 2 ~ 0.7. This 3D structure viewed in projection has no sig- 
nificant impact on recent mass estimates of A901b or the SW group components SWa 
and SWb. 
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1 INTRODUCTION 

Our current understanding of cosmology describes an intri- 
cate web of dark matter spanning the Universe dictating 
when and where galaxies form. These large-scale structures 
of dark matter have long been simulated in large N-body 
computations (Springel 2005), but due to the nature of dark 
matter, linking this feature of dark matter theory directly 
to observations has been a challenging task. 

Historically, the presence of dark matter as an essen- 
tial ingredient of the standard model of cosmology was de- 
duced indirectly from the rotation of galaxies or the sta- 
bility of galaxy clusters as the largest known cosmologi- 



cal objects. Nowadays, the dark matter density is inferred 
from the break-scale in the galaxy power spectrum, and 
by combining fluctuations in the cosmic microwave back- 
ground, the luminosity-density relation of supernovae Type 
la, the Lyman-a forest, or the bulk motions of galaxies. For 
an overview, see for example Dodelson (2003) and Peacock 
(1999) and references therein. Nevertheless, a direct detec- 
tion of dark matter particles as such is still an open issue. 

Weak gravitational lensing offers a unique way to chart 
dark matter structures in the Universe as it is sensitive to 
all forms of matter. Weak lensing has been used to map 
the dark matter in galaxy clusters (see for example Clowe 
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et al. 2006) with high resolution reconstructions recovered 
for the most massive strong lensing clusters (see for example 
Bradac et al. 2006). Several lensing studies have also mapped 
the projected surface mass density over degree scale fields 
(Gavazzi & Soucail 2007; Schirmer et al. 2007; Kubo et al. 
2009) to identify shear selected groups and clusters. The 
minimum mass scale that can be identified is limited only 
by the intrinsic ellipticity noise in the lensing analysis and 
projection effects. Using a higher number density of galax- 
ies in the shear measurement reduces this noise, and for this 
reason the Hubble Space Telescope (HST) is truly unique for 
this area of research. Compared to ground-based observato- 
ries, observations from space more than triple the number 
density of galaxies that can be resolved with a significantly 
enhanced imaging quality, permitting high resolution recon- 
structions of dark matter in the field (Massey et al. 2007; 
Heymans et al. 2008) and the study of lenses at higher red- 
shift. 

The lensing distortion experienced by a source galaxy 
depends on the total projected surface mass density along 
the line of sight (l.o.s.) and a geometrical factor that in- 
creases with source distance. The highest redshift galaxies 
therefore experience, on average, the strongest shear distor- 
tions. This redshift dependence can be used to recover the 
full three-dimensional gravitational potential of the matter 
density (Taylor 2001; Bacon & Taylor 2003) and has been 
applied to the COMBO-17 survey in Taylor et al. (2004) and 
the COSMOS survey in Massey et al. (2007). All these 3D 
mass reconstruction methods however require the use of a 
prior based on the expected mean growth of matter density 
fluctuations. Without the inclusion of such a prior, Hu & 
Keeton (2002) show that a Wiener filter approach is unable 
to reasonably constrain the radial matter distribution, even 
for densely sampled space-based quality lensing data. 

An advancement of the method of Hu & Keeton (2002) 
is presented in Simon et al. (2009), hereafter STH09, which 
details a non-parametric reconstruction method that di- 
rectly reconstructs the 3D matter density distribution. This 
is in contrast to previous methods which have mapped the 
gravitational potential. In broad terms one can view the 
gravitational potential as a heavily smoothed realisation of 
the density field. As an example, a density map of a single 
point mass will have a broad corresponding gravitational po- 
tential map that radially decays as 1/r in 3D (as In r in a 2D 
on the sky projection) from the position of the source. Di- 
rectly reconstructing the density field allows one to see the 
sources that produce the smoother gravitational potential 
through Poisson's equation and therefore this is ultimately 
what we wish to do. Historically, the 3D potential field has 
been mapped because the signal-to-noise was higher. Apply- 
ing a Laplacian to a noisy potential field appeared unfeasible 
on cluster scales. However, the work of STH09 showed that 
the Wiener filter method could also be applied in this regime 
to recover the matter-density field. Since the matter-density 
is the quantity we are interested in, and because the mass 
sheet degeneracy leads to unconstrained quadratic terms in 
the potential field, we now focus purely on the reconstruc- 
tion of the 3D matter-density field. Since this is the first 
application of this method, we shall pay close attention to 
the changes in signal-to-noise. 

Another new, computionally faster singular value de- 
composition technique has recently been presented in Van- 



derPlas et al. (2011). This paper proposes to use Eigenmodes 
in the reconstructed 3D mass density field, rejecting modes 
above a given noise level, thereby achieving a radial smooth- 
ing without a Wiener filter. On a deeper level, this is what a 
Wiener filter is doing as well, as it suppresses density modes 
according to their signal-to-noise, for which the prior is used 
to determine the expected signal in the reconstruction. 

The Wiener fiher employed by STH09 uses a ACDM 
prior to optimally smooth and reconstruct the 3D map. The 
smoothing does, however, have an unwanted impact as it 
(1) spreads out structures in radial (redshift) direction, (2) 
shifts on average the centre of structures to different red- 
shifts (here-after referred to as z-shift bias) and (3) scales the 
overall amplitude of matter density fiuctuations. We note 
here that while the radial spread and scaling of the ampli- 
tude of recovered perturbations was noted and detected in 
the 3D potential mapping approach, shifts in the position 
of structure were first pointed out in STH09 and therefore 
not accounted for in these studies. The information about 
how the true matter distribution is smoothed in the result- 
ing map is clearly defined by a radial point spread function 
(p.s.f). This radial p.s.f needs to be taken into full account 
for any correct interpretation of the map. In addition, as also 
outlined in STH09, the redshift resolution of all 3D lensing 
reconstruction will be low owing to the weak response of 
lensing to variations in lens redshift encoded in the broad 
lensing kernel. 

The STH09 method has been optimised for observa- 
tional data as it utilises the lensing information from all 
galaxies, whether they are with or without a photometric 
redshift estimate. It also has inbuilt tests for the effects of 
systematic errors in the lensing measurement, using an E/B 
mode decomposition. Crucially it also produces a 3D noise 
map so the significance of structures in the resulting maps 
can be assessed. 

In this paper, we present the first application of the 
3D matter density non-parametric reconstruction method, 
analysing the HST STAGES survey of the AbeU 901/902 
supercluster field, hereafter A901/902. We choose this field 
to test this method as it contains four known massive struc- 
tures at 2 ~ 0.165 detected at high significance in a 2D map 
(Heymans et al. 2008) and a fifth known structure in projec- 
tion at 2 ~ 0.46 (Taylor et al. 2004). In addition, STAGES 
is one of the largest HST surveys with one of the highest 
quality lensing data sets for this type of analysis. 

The 3D mass map is presented as a raw Wiener 3D 
reconstruction that reveals the expected low redshift reso- 
lution, based on our knowledge of the radial p.s.f and lens- 
ing kernel. Assuming the local maxima along each l.o.s. are 
a good estimate of the position of the peak density, we 
present 3D visualisations of the dark matter density in this 
field. Highlighting the most likely positions of structures in 
this way produces remarkable images that can be compared 
to other published 3D lensing reconstructions. However we 
voice strong caution in using any of these visualisations for 
scientific purposes as they present a very noisy estimate of 
the density field and contain no error information. More- 
over, structure along the same l.o.s. can overlap and thereby 
merge in the Wiener map to one single structure at a com- 
bined average redshift. 

For scientific purposes, we propose the use of raw 
Wiener reconstructions. We identify interesting l.o.s. across 
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the field from a 3D cross-correlation analysis of mass and 
light and then fully take into account the radial biasing 
and smoothing of the Wiener reconstruction to determine 
each lens redshift probability distribution. We then com- 
pare this probability distribution to over-densities of stellar 
mass along the same line l.o.s. and discuss our findings for 
the 3D matter distribution in the STAGES field. 

The paper is set out as follows. In Sect. 2 we review the 
STAGES data set and the 3D reconstruction method and 
present our results in Sect. 3. We discuss our findings and 
conclude in Sect. 4. As this analysis places new requirements 
on the accuracy of the weak lensing observations, the Ap- 
pendix details our tests of the observational analysis that 
ensure the results presented in this paper are robust. 



2 3D MASS RECONSTRUCTION OF STAGES 

In this section, we focus on the application of the STH09 
non-parametric 3D matter density reconstruction method 
to the STAGES A901/902 supercluster field, referring the 
reader to STH09 for details on the full theory that under- 
pins this method. This field has been the subject of sev- 
eral weak lensing analyses (Gray et al. 2002; Taylor et al. 
2004; Schirmer et al. 2007; Heymans et al. 2008), investigat- 
ing different aspects of the dark matter environment of this 
dense field. The most recent analyses have focused on data 
from the Space Telescope A90 1/902 Galaxy Evolution Sur- 
vey (STAGES, Gray et al. 2009). In this analysis, weak lens- 
ing distortions have been measured for over 60, 000 galaxies 
(~ 65 galaxies per square arcmin; Heymans et al. 2008) of 
which ~ 15% of the sample have accurate photometric red- 
shift and stellar mass estimates from the COMBO-17 survey 
(Wolf et al. 2004). Lensed galaxies without photometric red- 
shifts estimates are incorporated using magnitude bins for 
which the redshift distribution is estimated using Schrab- 
back et al. (2007). 

The 3D density reconstruction presented in this paper is 
the most complex lensing analysis of these data and places 
new demands on the accuracy and reliability of the weak 
lensing shear measurement and the accuracy of the photo- 
metric redshifts. In the Appendix, we therefore document 
the steps taken to confirm the reliability of the weak lensing 
shear catalogue that forms the basis of the 3D reconstruc- 
tion. 

At the very end of this section, in subsection 2.9, we 
give an executive summary of the main results to highlight 
the new results of this section. 



2.1 3D mapping in a nutshell 

In weak gravitational lensing, gravity-induced shape distor- 
tions of galaxy images are linearly related to the projected 
matter distribution along a line of sight (l.o.s.) direction 6. 
Let us start with a simple case in which we bin the complex 
ellipticities of all sources in a full shear catalogue onto the 
same regular grid. The grid covers the field-of-view of the 
survey, the mean ellipticity of the ith grid cell is denoted by 
e{0i); Oi is the l.o.s. direction of the ith cell. We arrange the 
binned ellipticities as one vector 



In the working regime of weak gravitational lensing, the 
convergence k{6) of light bundles is the projection of the 3D 
matter density fiuctuations S{r±,x), where rj_ is a comoving 
distance vector perpendicular to a fiducial l.o.s., and x is the 
comoving radial distance (Bartelmann & Schneider 2001): 
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Here c/Hq is the Hubble radius, /k(x) the angular diameter 
distance of x, f^m the mass density parameter and a(x) the 
scale factor at distance x- The function W{x) expresses the 
lensing kernel averaged for a distribution of sources in x: 
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Px(x) is the probability density of source distances. 

As Ansatz, we chop the matter distribution 5 into a 
set of A^ip matter slices or lens planes. Within each matter 
slice, 5 is assumed to be constant along a l.o.s. 6, so that 
the jth matter slice is represented by S^-'\6i), the average S 
in direction 0i (lens plane), and the width Axj ~ Xi+i ~ Xj 
of the slice; Xj defines the boundaries of adjacent slices. By 



means of 
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we arrange the grid values of 5 of every slice as a vector. 
Note that we assume a fiat sky where every lens plane is a 
tangential plane to the celestial sphere in the direction of 
the fiducial l.o.s.. 

Likewise, we arrange the k,{0) into a vector: 



K= (k(6Ii),k(6>2), 



(5) 



With this Ansatz, Eq. 2 can now be written as linear pro- 
jection of the (5-slices into the k grid: 



where 
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In the last step of Eq. 6, we denoted for convenience the 
weighed sum of d^'-''s by QS. 

In order to relate the convergence on the grid to the 
binned ellipticities of our survey, we have to transform con- 
vergence to shear 



7 = (7(0i),7(02),...)* 



(8) 



in the next step. On a fiat sky, the transformation is ex- 
pressed by the well-known Kaiser & Squires (1993) relation, 
which is written down for a regular grid here (Hu & Keeton 
2002): 



7 = Pj^Qd 



(9) 



where we have for the convergence-to-shear transformation 
matrix 



(e(0l),e(02). 
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with 9ij = Oi~ Oj and A being the solid angle of a grid pixel. 
Note that for the 2D vector we are using the commonly 
employed complex notation. 

So far, nothing has been said about the noise, 

n= (n(0i),n(02),...)* , (H) 

which is necessarily included in e owing to the intrinsic shape 
of the sources and the non-uniform distribution of sources 
resulting in occasional cases of grid cells with zero sources 
(for example in a masked region): 

e = f + n = P-,^QS + n. (12) 

Of n only the statistical properties are known: it is assumed 
to be uncorrelated to S, has a vanishing statistical aver- 
age and the covariance N — inn'). In the case of missing 
sources, the noise n{Oi) associated with the zth grid cell is 
infinite. 

In 3D lensing, we have a survey that can be split into 
various subsamples of sources with distinct and known red- 
shift distributions. Now each source subsample provides a 
different grid e^, which every time is a projection of the 
same matter cube 8 but with a different set of weights 
and noise grids rii. Therefore, employing a compact nota- 
tion, we now have 

e = (ei,e2,...) 

= P^«(Qi<5,Q2<5,...)-f (ni,n2,...) 

= P^„Q(5 + n . (13) 

The map making algorithm utilises this relation to de- 
fine an estimator for 5 given e and the noise properties N. 
Although Eq. 13 has no exact solution due to the unknown 
noise component n and the non-invertible P^k (mass sheet 
degeneracy), one can find an optimised solution, i.e. a linear 
filter H of e, that minimises the average deviation between 
the estimated minimum variance S-^v = He and the true S, 
i.e. 

(II .5n.v -<5 11") =min. (14) 

(By II X |p= x^x we mean the Euclidean norm of x.) 

The solution to this optimisation problem is in our case 
(e.g. Zaroubi et al. 1995): 

H = [al-f SR'''N"^r]~^SR'''N"\ (15) 

where R = P7KQ, and a is an additional constant, usually 
between a G [0, 1], that allows to change the signal-to-noise 
in the resulting map by regulating the smoothing of the fil- 
ter (Saskatoon filter; Tegmark 1997). The signal covariance 
S = (ee^) — N has to be given to the filter as prior and is 
computed for a particular fiducial cosmological model. The 
practical challenge is the computation of 5mv = He within 
a reasonable time, see the Appendix of STH09 for details. 

To evade the problem of the mass-sheet degeneracy, we 
restrict ourselves to solutions for which the average of (5mv 
vanishes for every redshift slice individually. 

2.2 Reconstruction prior and grid parameters 

The prior can be implemented in the form of a radial or a 
transverse filter; a transverse filter varies as a function of 
angular wave number i (STH09). For S, we choose a radial 
Wiener filter 



(5«(0,)5('='(0O) = / d«P]'^(^)|F(e.^)|Vo(^0) ,(16) 
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where F{x) — 2x~^ Ji{x) is the Fourier transform of a top- 
hat window, Jn{x) a Bessel function of first kind and Pg''\() 
is the power of matter fluctuations, projected along the 
width of the ith matter slice. We approximate our square 
grid pixels by circular pixels of identical area A = ttOj. 

As a fiducial matter power spectrum we use the approx- 
imation of Smith et al. (2003) and Limber's equation for con- 
verting 3D to 2D power on the sky. The fiducial cosmology 
is a ACDM model (adiabatic fiuctuations) with Q.m = 0.24, 
rib = 0.0416, r2A = l-r2m and Ho = /i 100 km s~^ Mpc~^ 
with h = 0.732. The normalisation of the matter fluctuations 
within a sphere of radius 8ft~^Mpc at redshift zero is as- 
sumed to be og, — 0.76. For the spectral index of the primor- 
dial matter power spectrum we use jig = 0.96. These values 
are consistent with the third-year WMAP results (Spergel 
et al. 2007). 

The smoothing of the Wiener fllter can be regulated by 
rescaling the noise covariance of the data via a: For a = 
Wiener filtering is switched off, i.e. the prior S cancels out 
at the expense of the lowest signal-to-noise of the recon- 
struction but for the benefit of an unbiased estimate of S. 
For a = 1 the Wiener fllter boosts the signal-to-noise to a 
maximum at the expense of heavy smoothing in the radial 
direction along with a strong z-shift bias and loss in redshift 
resolution. We find a = 10"'^ provides a good balance to 
maximise signal-to-noise and minimise bias, which we fully 
account for below. Note, a radial prior does not produce 
transverse smoothing so an additional transverse smooth- 
ing is required for every lens plane. For post-reconstruction 
smoothing we use a Gaussian kernel of arms ~ 40 arcsec, 
which is the size of several grid pixels. 

The number of matter lens planes used here is A^ip = 22, 
covering a range from 2 = ... 1.15 with an associated slice 
width of Az — 0.05. There is one exception for the closest 
plane which has Az = 0.1. The width of the slices is finer 
than the radial resolution we can possibly hope for. Never- 
theless, the large number of lens planes was chosen to avoid 
numerical instabilities in the reconstruction. The lens planes 
are gridded with an angular resolution of 128 x 128pixel^, 
one pixel has a square size of 16.34 arcsec. The vector S 
hence has 128^ x A^ip « 3.6 x 10^ de grees of freedom. 

2.3 The radial point spread function 

As outlined in STH09, the reconstruction algorithm does not 
provide an unbiased representation of the true radial matter 
distribution. The bias can however be quantified. 

Ignoring the noise in a particular reconstruction, on av- 
erage the map Smv will be a smoothed image of the true 
matter distribution <5: 

(<5mv> = H(e> = HP.,«Q(<5 + n> = HPy^QS , (17) 

where (...) is the ensemble average over all source noise real- 
isations while keeping the actual matter distribution S con- 
stant. Hence the matrix B = HP^y^Q holds all information 
on the smoothing involved in the reconstruction process. In 
the optimal case, one would have B = 1. Here, we do not 
have B = 1, though, but still B contains all information 
needed for a quantitative assessment of the resulting map. 
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Unfortunately, as discussed in STH09, simply applying 
the inverse of B to 5mv, even if B was regular, would undo 
the Wiener filter, leaving behind a heavily oscillating recon- 
struction with a significant noise level because low signal- 
to-noise modes in the map would no longer be suppressed. 

Here, the full matrix B has of the order of 10^ x 10^ = 
10^° elements. It turns out, however, that in practise most 
of its elements are negligible because the devised Wiener 
filter smoothes the matter distribution only in the radial 
direction. 

Concretely, consider that we have just one point source 
overdensity inside the entire 3D density field S, represented 
by one unity-valued grid pixel in direction Oi on the jth 
lens plane. To lift the mass-sheet degeneracy, the algorithm 
assumes that the density contrast averaged over every lens 
plane vanishes. Hence, to stay inside the space of considered 
<5-configurations all other pixels of the jth lens plane have 
to be set to -(127 x 128)"^ for our (128pixel)^ grids; pix- 
els of all other lens planes i ^ j are zero. We denote this 
special density contrast configuration by e"/^ ■ Any arbitrary 
other configuration <5, with zero average contrast for each 
lens plane, can be represented as weighed sum of ep''s as 



■ Cijef ^ with coefficients Cij. The reconstruction al- 
gorithm is linear so that the resulting noise averaged Wiener 
map is a linear combination of all Be'"*', i.e. 



H(P.,«Q<5 + n) ^G5 + H{n) C^Be, 



(18) 



the vector n denotes the shear noise in a particular realisa- 
tion. 

Essentially the original matter point source e^p is 
smeared out in the l.o.s. direction Oi, but not in the trans- 
verse. In STAGES the radial smearing will be roughly un- 
changed when varying 9i inside the same lens plane i. We 
checked this by computing Be^' for a set of A'los = 10 x 10 
l.o.s. directions for all lens planes covering the inner 80% of 
the field. Therefore, instead of keeping the full Be^' for all 
ep', we store only the average l.o.s. profile of Be'"*^ along 6i 
and average over several Oc 



P 
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where [A, Oi] denotes the linear operation that takes only 
elements from the 3D grid A that are along the l.o.s. direc- 
tion Oi and arranges them as an A^ip element vector in order 
of ascending lens plane redshifts. Note that Bep^ is calcu- 
lated in practise by applying the reconstruction algorithm 
to P^^Qe*^' so that in total iYi^, x Nip = 2200 individual 
reconstructions were performed. 

Fig. 1 shows the resulting radial p.s.f. for a selected 
set of lens planes. Radially located peaks in the true matter 
distribution will have a cigar-like appearance in the map due 
to the evident width of the p.s.f. The spread will be reduced 
for better signal-to-noise in the data, but will always be 
present to some degree when a Wiener filter is employed. 

2.4 ^-shift bias 

An obvious application of the Wiener map is to estimate 
the redshifts of 3D peaks by the local maximum of the cigar 
density profile or to compare the redshift estimate to other 
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redshift 

Figure 1. Average radial smearing for different lens plane red- 
shifts, shown as small numbers attached to the profiles. Note that 
the local maximum of a p.s.f. generally does not coincide with the 
lens plane redshift. 




0.4 0,6 0.. 

true redshift 

Figure 2. Z-shift bias of the radial filter used in the reconstruc- 
tion. Owing to the Wiener filtering, structures do not appear at 
the correct redshift in the lens planes as indicated by the data 
points. Using this figure the 2-shift bias in the reconstructed 3D 
map can be corrected for (see text). The error-bars denote the 
r.m.s. -variance of the bias over the field-of-view of STAGES. Open 
squares: bias in the density contrast map; filled points: bias in the 
S/N map. 



maps of the galaxy number density or luminosity. However, 
as can be seen from Fig. 1, generally the maximum of the 
p.s.f. does not coincide with the true redshift of structures, 
which we hereafter refer to as z-shift bias. To determine the 
extent of z-shift bias we use Fig. 1 to determine where, on av- 
erage a point source mass is recovered in the reconstruction. 
The result is shown in Fig. 2, which displays the redshift 
of a p.s.f. local maximum (apparent redshift) and error in 
relation to the original point source redshift in e'f\ We use 
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the data points in this figure to construct a new Wiener map 
free of the z-shift bias in the foliowing way: We take a set 
of new lens planes with "true redshifts" and fill for every 9i 
these planes with values from the uncorrected Wiener map 
at the "apparent redshift" given by Fig. 2. 



2.5 Quantifying noise and systematics 

To assess the S/N of features inside the map, we divide map 
values by the variance observed in 10^ map noise realisations 
before z-shift bias correction. Noise maps are realised by 
randomly rotating the orientations of sources and rerunning 
the reconstruction algorithm. During this process, we also 
average the noise covariance N of density pixels along the 
same l.o.s.. 

Dividing by the noise variance inside a lens plane 
changes the shape of the p.s.f. and the z-shift bias, see open 
squares in Fig. 2. In the case of STAGES, the S/N-map 
shifts all p.s.f. peaks essentially to the middle of the map 
(in redshift), which can be explained by larger levels of noise 
variance at low and high redshifts. Due to the strong impact 
of the z-shift bias, the S/N-map is in essence thus only a 2D 
map stretched out over the entire redshift range. For that 
reason, we will use it in the following only as a 2D projection 
to quote the maximum S/N expected in a l.o.s. direction. 

To test the impact of systematic errors in the lensing 
shear measurement, B-mode maps are produced. These are 
obtained by rotating the source ellipticities by 45-degrees 
followed by a repetition of the reconstruction process. Ide- 
ally, these maps should be pure noise maps as a gravita- 
tional tidal field is unable to generate a B-mode distortion 
by shearing galaxy images. Instead, they arise from system- 
atics related to the source shape p.s.f. correction procedure, 
the intrinsic clustering of sources (Schneider et al. 2002) and 
intrinsic galaxy alignment correlations of physically close 
sources or shape-shear-correlations (Heymans et al. 2006, 
and references therein). Moreover, as discussed in STH09, 
matter fluctuations outside the field-of-view but close to 
the field-of-view may be capable of generating large-scale 
B-modes within the reconstruction volume. Therefore, an 
entirely vanishing level of B-modes is never guaranteed per 
se. To partly suppress this effect, we take a reconstruction 
field-of-view somewhat larger than the actual coverage of 
the survey and later remove the area not covered by obser- 
vations. This allows unknown sources of shear to be outside 
the field-of-view and to be properly accounted for in the 
model reconstruction. 

The peak and trough statistics in Sect. A3 shows 
that the B-mode map is consistent with noise realisations, 
whereas the E-mode map is inconsistent with noise at high 
confidence as it should be. This leads us to the conclusion 
that any systematic errors are low and within our statistical 
errors. 



2.6 Redshift uncertainties 

The main drawback of the Wiener method is that the radial 
spread of the p.s.f. hampers the ability to attach a redshift 
to a cigar-shaped profile in the map. This would be true 
even if a p.s.f. was very narrowly peaked as the uncertainty 
of the estimate still needs to be determined by accounting 



for the noise in the map, which is not directly visible in the 
Wiener estimate for the 3D matter density contrast. 

What really determines the redshift accuracy is the abil- 
ity to distinguish a radial profile p*'^ from p'-*', Eq. 19, de- 
spite of the noise superimposed onto the average profile in 
a concrete realisation. 

To get a handle on the principle redshift accuracy, let 
us consider a simple case where the map density profile in 
the l.o.s. direction is dominated by one single structure on 
the ith lens plane, i.e. <5'''' — (5oP*'^ + n; So is the grid pixel 
density contrast and n the noise vector in that particular 
realisation. To pin down the redshift of the structure from 
the noisy vector d''-* means whether we can statistically dis- 
tinguish, given the data, a model (a) in which the mass is 
presumed to be on lens plane j from the correct model (b) 
where the mass is on lens plane i. Accepting model (a), we 
would estimate the density S'o — XSo at redshift Zj by min- 
imising 



N 



(20) 



We are assuming a Gaussian noise model with covariance 
N — {nn'^) here. Since the source noise is roughly Gaussian 
and the Wiener filter linear this is a valid assumption here. 
The minimum solution for the previous equation is 



(21) 



so that the oi the best-fit of model (a) has 

X^J = (Amin^op'^' - ^ ( A,„i„<Sop''' - , (22) 

which on average for all realisations n is 

ixl) = Ax?,+iV,p, (23) 



(25) 



Therefore, to distinguish model (a) and model (b), we have 
to decide when {xij) is inconsistent with a x'^-distribution 
with A''ip degrees of freedom, i.e. the x'^-distribution of model 
(b). Since model (b) has Amin = 1 and (Xij) ~ Nip, the 
difference in (x^) is simply the above Ax?j. 

To get fiducial values of Ax?j for STAGES, we model 
So as the centre of a singular isothermal sphere (SIS) with 
mass profile (e.g. Bartelmann et al. 2001) 



Ssis{d) = 



4tv c^a(Xh) 

3 H^n m 



1 



(26) 



/k(xh)e ' 

where Xh is the SIS comoving distance and crj the velocity 
dispersion of the SIS. The SIS is smoothed it about the 
centre with a Gaussian kernel of size o"rms = 40 arcsec to the 
Wiener map transverse smoothing scale: 



Soixh) 



(2^)3/2 c2a(xh) 



J J Q i im 



1 



(27) 



AXh/k(Xh)ci"rms 

Note that we also have to divide the smoothed Ssis (SIS den- 
sity contrast integrated along the l.o.s.) by the width Axh of 
the matter slice it resides in, to comply with the definition 
of the mean density contrast in our mapping algorithm, Eq. 
6. 
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Figure 3. Expected uncertainties of redshifts of structures with 
a SIS-like density profile as function of M200 mass (units of 
lO^^M© /i-i). The thick solid black Une shows the average red- 
shift (true redshift: x-axis; apparent redshift: j/-axis) of a SIS (no 
z-shift bias). The upper and lower lines bracket the 68% con- 
fidence region for the scatter of redshifts inside the map. The 
width of the confidence interval depends on the SIS mass and its 
true redshift. 



Based on the Ansatz outlined above, Fig. 3 displays the 
expected accuracy with which redshifts of structures can be 
determined from a Wiener filtered 3D map. In this plot, 
masses are defined as SIS-M200 masses, i.e. the mass con- 
tained inside a radius r-2oo within which the mean density 
equals 200 times the critical density: 



10GH{z) 



(28) 



6.58 X 1O"M0/i"^ 



103 kms- 



Ho 



For comparison, the two most massive clusters in the field, 
A901a and A901b, have approximately a mass of M200 ~ 
2 X lO"M0/i"^ (Heymans et al. 2008). The figure adopts 
the average radial p.s.f. of STAGES. For a STAGES-like 
survey we find the redshift of a massive galaxy cluster of 
M200 ~ 10^^ M0 could still be measured with accuracy 
of A.Z f» ±0.08 (±0.15) at 2 ~ 0.6 (0.8), whereas a cluster of 
mass with M200 ~ 3 x lO^^M© /i"^ has Az « 0.05 at z ~ 0.2 
but already Az « 0.15 at z ^ 0.6. For masses smaller or 
equal than M200 ~ , the accuracy quickly de- 

teriorates with redshift, at most imposing upper limits on 
the redshift. This lack of 3D resolution for lower mass ha- 
los poses the natural question of how useful a 3D lensing 
reconstruction is, which we discuss in Section 4. 




Figure 4. The minimum M200 mass required to have a detection 
of a SIS located at a certain redshift (x-axis) above a certain 
significance level (different lines) anywhere in the Wiener map of 
STAGES. 



[So{Zt)p' 



jG{l,...,Mp} 



(29) 



along the SIS l.o.s. in the Wiener map; STAGES survey 
parameters are assumed. For a given Zi, M200 is changed 
until Ran reaches a desired S/N limit. This limit is easily 
found due to the scaling A'hoo oc Rln, which follows from 
Eqs. 28 and 27. 

For A901a/b at 2 = 0.165 we expect a S/N between 
6 — 7(7 in the Wiener map, which indeed is found. We can 
also see from this plot how well suited the STAGE field is for 
this type of analysis as we are primarily sensitive to struc- 
ture at low and moderate redshift for M200 ^ IO^^Mq/i"^. 
To reconstruct mass at higher redshift requires a significant 
structure to lie in the field. 



2.8 Product Map Mass and Stellar Mass 
comparison method 

For a further analysis, we will compare the lensing mass map 
to the distribution of stellar mass inside the data cube. The 
stellar mass map is derived directly from the COMBO-17 
galaxy catalogue, which is binned onto the 3D grid used for 
the mass reconstruction. We compare the two 3D gridded 
distributions on a pixel-by-pixel basis with the aim to find 
local matches and discrepancies between the light and mass 
distribution. 

For the mass- light comparison, we restrict the light map 
to a selected redshift range: 



2.7 Detection limit for mass concentrations 

In order to gauge the mass detection limit of the STAGES 
Wiener map, we have produced Fig. 4. For this plot, we 
place the central pixel So{zi) of a SIS at lens plane redshift 
Zi and determine the maximum S/N ratio. 



(30) 



where the indices l,k = 1 . . . Nip, I ^ k choose the redshift 
range in question; Bi denotes a vector with Nip elements 
being zero everywhere except for the ith element. The range 
is 2 £ [zk, zi], where Zi is the redshift of the ith lens plane. 
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We understand the lensing map as a convolved under- 
lying matter distribution, 



= PS 



(31) 



The convolution of the mass map is an unavoidable by- 
product of the Wiener filter reconstruction. Therefore, in 
contrast to the light map, for the mass map we do not have 
the choice of selecting a certain redshift range. We there- 
fore cross-correlate stellar mass with the full lensing mass 
map. This is reasonable because distinct lens planes are well 
separated enough in physical distance to have no relevant 
astrophysical correlation, i.e. for i ^ j 







(32) 



Our analysis tool is a product map where the product 
is taken of the stellar mass and total mass density profiles 
in radial direction: 



Vlp I 



which on average is 



(33) 



(34) 



where we define the weights Wi = [p''']*ei. From the latter 
it becomes obvious that Eq. 33 can be made an estimator 
for the weighted l.o.s. product of mass and light distribu- 
tion within the selected redshift range [zfc,^;], namely by 
the product map (PM) 



PM(0,fc,O^^^^^; 



(35) 



which for statistically independent lens planes has the aver- 
age 



PM{e,k,l) 



'S' 



(36) 



We use product maps to obtain candidates for l.o.s. direc- 
tions for which features in the mass map could be related to 
stellar mass (PM > 0), directions that are significant in the 
mass map but apparently uncorrelated to the stellar mass 
map (PM = 0) or even anti-correlated (PM < 0). The can- 
didates will be subject to further modelling of the radial 
density profile to determine the probability distribution of 
mass along the l.o.s.. Sect. 3.2. 

The product method has a limited radial resolution due 
to the broad radial p.s.f. p'^'K For illustration, consider a 
l.o.s. that has exactly one non-zero pixel Ai in Si on the ith 
lens plane. On the other hand, the matter distribution may 
have the only non-zero pixel Am on the jth lens plane, hence 
at a different redshift. Despite an obvious mismatch between 
light and mass, this case will contribute to the product map 
with AiAmieilV''' (Eq. 33), if p'-'' is broad enough to reach 
the ith lens plane. 

In Fig. 5, we show the response, or signal, in the prod- 
uct map to a stellar mass point source at discrete redshifts 



o o 
a- 6 





7033 ^^'^^ 

/ 0.58 \ 










0^ 










1.06 




.11.11 K mop 








— 3-D map 







0.2 0.4 0.6 0.8 1 

redshift mass map 

Figure 5. A matter density fluctuation in the mass map still 
contributes to the product map even if it is at a different redshift 
(x-axis) as a fluctuation in the stellar mass map (redshift of filled 
triangle). This is an effect of the broad radial smoothing kernel 
the Wiener map is subject to. In an ideal world, the solid lines 
would be strongly peaked about the triangles, i.e. [p'')]'ej 5^. 
Thicker lines belong to higher stellar mass redshifts. For compar- 
ison, the dashcd-dotted line shows the response we would expect 
for a purely 2D (convergence) map; the response curve is inde- 
pendent of the stellar mass redshift in this case. 



z = 0.05, 0.33, ...1.08 (shown with triangles) from a lens- 
ing mass spread over the full redshift range < z < 1. Mass 
along the full redshift range will contribute at some level to 
the product map, even though it is not associated with the 
stellar mass point source. This shows that our interpretation 
of the product map is non-trivial as there is always consid- 
erable confusion present for the STAGES map. However it 
does show that we have relatively little confusion from low 
redshift lensing mass features to stellar mass redshift range is 
z > 0.5. This is important as it means one can, in principle, 
distinguish in the product map the supercluster structure 
at 2 ~ 0.165 from structures in the higher redshift regime. 
The response could be further improved by using less ra- 
dial smoothing in the Wiener map by means of devising a 
smaller tuning parameter. This, however, also decreases the 
signal-to-noise in the product map. 

A cross-correlation of a stellar map with a 2D conver- 
gence map, on the other hand, always produces the same 
response irrespective of the stellar mass redshift. This re- 
sponse is shown in Fig. 5 by the dotted-dashed blue (am- 
plitude was rescaled for plotting purposes). Although the 
cross-correlation with the 3D Wiener map displays some im- 
provement over a cross-correlation with a convergence map, 
it is unclear how the still broad response can be properly 
dealt with. For instance, a massive mass peak at low red- 
shift could generate a strong response with a large stellar 
mass peak at high redshift. 

To advance the 3D product map, it is conceivable to de- 
fine a radial smoothing of the light map in order to optimise 
the response in the product map by substituting by an 
appropriate vector. In addition, one could properly account 
for the z-bias so that the maximum response is at the loca- 
tion of the stellar mass peak or trough, which is not exactly 
the case here (Fig. 5). We do not pursue such a strategy 
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for this work, though, as we are not expecting a significant 
improvement for STAGES, fn particular, we do not use a 
3D product map to draw strong conclusions about radial 
matches of stellar mass and matter mass. Instead, we solely 
employ the product map in combination with the Wiener 
S/N-map to identify interesting l.o.s. to be further analysed 
by a purely lensing method Sect 3.2. 

Note that smoothing the light map with the Wiener 
map p.s.f. before multiplication, i.e. set — p^^\ results 
in a worse response than the one adopted for this analysis; 
the response would be the product of pairs of Wiener p.s.f. 
equalling their overlap, i.e. Wi — [p'''']'p'''' in this case. 

Finally, the signal-to-noise of the product map is as- 
sessed by dividing Eq. 35 by the variance seen in noise reali- 
sations, which are obtained by randomising the positions of 
galaxies in the M, map. Note that the signal-to-noise levels 
in the product map can become considerably larger than in a 
pure lensing map because for uncorrelated noise of the lens- 
ing and Af, map the signal-to-noise in the individual maps 
just multiplies. 

2.9 Executive Summary 

We have reviewed the theory presented in STH09, then 
presented details of the extensions and new analysis tools. 
These tools were developed for the first application of this 
method to real data. Here we highlight the new results pre- 
sented above for different readers. 

For readers optimising future survey designs, the most 
relevant results are contained in Figs. 3 and 4, which show 
the best 3D mass resolution that can be attained on the con- 
sidered angular scales with space-based quality data with 
a galaxy number density of 65 galaxies per square arcmin 
matched by 10 galaxies per square arcmin with ground- 
based photometric redshift estimates with a redshift error of 
az/{l + z) = 0.02. Note these photometric redshifts are lim- 
ited at z < 1.4 in the absence of near infra-red photometry to 
extend reliably to higher redshifts. This survey characterisa- 
tion is similar to planned deep imaging areas of space-based 
surveys with matching photometric redshifts from shallower 
multi-colour optical surveys such as PanSTARRS, DES or 
VST-KIDS. These reveal rather relatively poor expectations 
for attaining mass resolution in 3D using lensing alone for 
anything other than the most massive structures (at least 
several lO^^AI^h'^). It therefore motivates future surveys to 
invest in sufficiently deep photometric ground-based surveys 
such as LSST and matching deep near infra-red photometry. 

For readers focused on applying the STH09 method to 
a new data set. Sect. 2.3 and 2.4 detail how to make the non- 
trivial measure of Z-shift bias and the radial p.s.f from real 
data that was not discussed in STH09. Fig. 2 also shows how 
the application of standard 2D mass mapping noise estimate 
techniques to a 3D mass map biases the redshift measure 
with the implication that a 3D S/N map, which has not 
previously been calculated in the literature, is essentially 
a 2D map for lensing surveys like STAGES. Section 2.8 
therefore outlines a new product map method that enables 
one to compare significant peaks detected in the essentially 
2D S /N mass to the 3D stellar mass distribution and the 3D 
lensing density distribution. 

We show that the product map offers theoretical im- 
provements over a cross-correlation of a 3D stellar mass 



map with a 2D convergence map, but it remains unclear 
at this point how to clearly exploit this map. Nevertheless, 
we employ the product map to select interesting l.o.s. to be 
analysed further in the following. 



3 RESULTS 

In this section, we present the results of our 3D recon- 
struction of the STAGES A901/902 supercluster field and 
our analysis of that reconstruction. We start by presenting 
the raw Wiener reconstruction in Fig. 6, splitting the field 
into two parts: the northern half of the field, comprising 
A901a/b (left column), and the southern part, comprising 
A902 and the SW group (right column). In the upper pan- 
els we present the 3D density map projected on the plane 
of the sky showing with high significance the four structures 
in the supercluster. This projection shows the same result 
as the 2D reconstruction in Heymans et al. (2008) and we 
overlay signal-to-noise contours obtained by comparing the 
mass map to noise realisations. Dividing by the noise in the 
3D mass map, changes the z-shift bias in the new resulting 
map, rendering it essentially a 2D map for STAGES. For 
this reason, the S/N map is shown only in 2D as maximum 
along a l.o.s. direction. In the lower panels, we present the 
3D map now projected along the declination. This recon- 
struction shows the redshift-density distribution convolved 
with the the radial p.s.f. of the map, shown in Fig. 1, cor- 
rected for z-shift bias, shown in Fig. 2. As this radial p.s.f is 
so broad, structures are radially spread out in the map and 
overlap in redshift. Nevertheless, in cases where structures 
do not lie in close projection, once the z-shift bias is cor- 
rected local maxima are unbiased estimators of the redshift 
of mass concentrations. Whilst unbiased however they have 
limited redshift accuracy, refiected by Fig. 3. 

In order to better visualise the 3D matter distribution 
from the raw Wiener map, we make the general assumption 
that the local maxima are good estimators of the peaks in 
the density field. We then create a 3D map retaining only 
the density pixels that are a local maxima along each l.o.s.. 
We then smooth the visualisation with a Gaussian kernel 
(0.2 arcmin)^ in r.m.s. size on the sky and 0.025 in red- 
shift with the result shown in Fig. 7. This cleaning process 
produces a striking image of the 3D density field revealing 
structure at high redshift (marked SWc and b/3) behind 
the foreground supercluster (marked A901a, A901b, A902, 
SWa/b). It is this realisation of the noisy 3D density field 
that should be compared with other 3D mass reconstructions 
in the literature (Massey et al. 2007; Taylor et al. 2004). A 
one-to-one comparison is still not easy, though, since other 
reconstructions show the estimated gravitational potential. 
We strongly caution the reader over the direct use of 3D 
lensing maps, however, as they are noisy estimators of the 
redshift distribution and contain no error information. In 
particular, owing to the broad smoothing kernels, structure 
along the same l.o.s. can merge in this visualisation to pro- 
duce a density profile with a local maximum at an intermedi- 
ate redshift. As our detailed analysis of the raw Wiener map 
will show however, after the radial p.s.f has been correctly 
accounted for, the higher redshift structures picked out in 
this visualisation are still seen with a high probability. 
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Figure 6. Different projections of the Wiener mass map (top row: projection onto the sky; bottom row: radial projection of the area in 
the top row; 2-shift bias corrected) The colour scales shows the maximum density contrast in the map along the projection axis. The 
contours in the top panels correspond the S/N levels 2.5(t, 3(t, 4(t, 5(t and 6a" as determined from noise realisations of the data. Contours 
in the bottom panels are density contrast levels included to highlight the structure in the radial map. Structures appear radially stretched 
due to the impact of the Wiener filter smoothing kernel. 



3.1 3D comparison of mass and light 

One of the main science drivers for 3D mass mapping is to 
investigate how well the total mass traces light or stellar 
mass in 3D. We use the product map technique outlined in 
Sect. 2.8 to compare total mass and stellar mass in redshift 
slices and identify interesting l.o.s. in our 3D map. 

We use 17- band COMBO-17 stellar mass M, estimates 
for the full galaxy sample down to i? < 24 mag. The M*-map 
is produced by summing the stellar masses Mi of galaxies 
within the same angular bin and lens plane (pixel). Bins 
containing no galaxies have a stellar mass of zero. This ma- 
trix of binned masses is divided by the average pixel mass, 
M = Mi/Np where A^p is the total number of unmasked 
pixels. All pixels of the same l.o.s. are flagged unmasked if 
there is at least one pixel along the l.o.s. containing one or 
more galaxies. The majority of masked pixels are found at 
the edges of the field. The effect of masking in the central 



part of the map is hence small. The normalised M,-map of 
the lens plane is in the unmasked regions, up to a constant, 
an estimator for the M* density contrast, M{0)/M — 1. In a 
next step, the average M{6) / M is subtracted from all pixels. 

In Fig. 8 we present the Total-Stellar Mass Product Map 
of the STAGES field within four redshift shces z £ [0, 0.25], 
2 G [0.25,0.50], z e [0.50,0.70] and z G [0.70,1.00], where 
the redshift bin width is significantly larger than the pho- 
tometric redshift error a, = 0.037. The colour-scale shows 
the product map level normalised by the variance in noise 
realisations; black through purple regions have negative val- 
ues and blue through red have positive values. We can see 
the majority of the density map is positively correlated with 
the stellar mass map. Overlaid in red are S/N contours of 
the density map. As shown in Sect. 2.4, the S/N maps are 
almost completely spread out in redshift and so we use the 
same S/N map in all redshift slices to identify regions in 
the plane of the sky where we find a significant detection of 
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Figure 7. 

3D Dark Matter Visualisation. Representation 
of the local l.o.s. maxima in the raw 3D Wiener 
map of STAGES, obtained from Fig. 6. This 
visualisation of the 3D density field is an noisy 
estimator of the true distribution of mass con- 
centration. The radial spread of structure does 
not reflect the redshift uncertainty. 
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Figure 8. Total-Stellar Mass Product maps obtained by multiplying a stellar mass map with the Wiener mass reconstruction within 
four rcdshift slices (top left for ^-rangc). Intensity levels in the background show the signal-to-noise in the product map with black 
regions having the most negative values and red the most positive levels. Signal-to-noise mass density contours are shown in red with 
levels 2.5cr, 3a", 4cr, 5(7 and 6cr, depicting the maximum significance along a l.o.s. in the lensing mass map. Note the source of the signal 
may be at any redshift as the S/N maps are almost completely spread out in redshift. White contours highlight regions of high stellar 
mass Mh. = 10^", 2 X 10^"^, 5 X 10^'^, IO^^Mq, in the M^, map smoothed with a 1.2 arcmin Gaussian kernel. The map pixel volume 
increases towards greater redshift. Open diamonds are reference points marking prominent bright cluster galaxies in A901a/b, A902 and 
the SW group. Numbers on the x- and y-axis denote the right ascension and declination, respectively (J2000.0). 
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mass. The 3D stellar mass distribution is shown with white 
contours. We find a clear correlation of mass and stellar mass 
in the z — 0.165 A901/902 supercluster in the lowest red- 
shift panel (upper left) and also the CBI cluster ai z = 0.46 
(upper right). 

The interpretation of this map is however non-trivial 
as the radial p.s.f. of the map causes overdensities in the 
mass map to contribute to the product map even when they 
are at a different redshift to the concentration in the stellar 
mass map (Sect. 2.8). We can therefore only use this prod- 
uct map to identify interesting l.o.s., which we then follow 
up with a detailed analysis that incorporates knowledge of 
the radial p.s.f. of the Wiener map. We identify interesting 
l.o.s. as follows. First, we only consider pixels in the total 
mass map that are detected with a S/N > 3 in projection. 
We then choose all pixels where we find a peak positive or 
negative correlation of significance greater than la with the 
stellar mass. The majority of the lines of sight above the 3cr 
threshold exhibit a positive correlation in the product map. 
According to the noise peak statistics in Appendix A3, how- 
ever, we expect at least some of the matches below < 3. So- 
to be confusions of random lensing map fluctuations with 
the 3D stellar mass map. We find that the cases of negative 
correlation in the lensing map with S/N > 3 are all behind 
foreground structures A901a, A901b, A902 and SW, which 
are likely to be produced by the leakage effect, outlined in 
Sect. 2.8. With the exception of A901a, these negative corre- 
lations have a product map signal-to-noise of < 3cr and are 
therefore statistically insignificant. They are all included in 
the analysis of the selected l.o.s. with positive correlations in 
the foreground, and therefore automatically included in the 
following analysis. Our ten chosen l.o.s. are listed in Table 

1 with identifiers based on the known foreground z = 0.165 
structures in this group, A901a, A901b, A902, SW, and also 
the higher redshift CBI cluster at z = 0.46.^ 

Before continuing with our detailed line of sight mod- 
elling of the ten structures marked in the product map, 
we briefly discuss other strong features. In the two high- 
est redshift slices, we find two stellar mass overdensities 
at M, = 1O"M0 (with 1.2 arcmin smoothing). In the 

2 £ [0.50, 0.70] bin, we see a strong positive correlation PI of 
total mass and stellar (shown white) where the total mass is 
detected at 2.5cr. In the z-bin z £ [0.70, 1.00], we find signif- 
icant anti-correlation P2 (~ 6a in the product map), which 
also is a stellar mass overdensity matched by a lensing mass 
underdensity. The M« overdensity region is found to be over- 
lapping with gaps produced by three stars in this area, which 
may be interfering with the lensing map. Moreover, as shown 
in Fig. 4, at these high redshifts we would only expect to re- 
liably detect haloes of mass M200 ~ 2 x IO^^Mq, a higher 
mass range than what we would expect for these stellar mass 
overdensities. The same argument holds for the feature P3. 
These stellar mass overdensities therefore demonstrate the 
effect of noise in our mass maps where we are detecting mass 
with lensing only at the 2.5a level. We also note that we find 



^ Note that SWa/b, galaxies belonging to a known group at 
z = 0.165, did not make it into the table because the lensing 
significance in direction of PM peaks corresponding to SWa/b is 
too low; for SWb it is 2.5(T. However, SWb lies within the range 
of the SWc pencil beam analysed later on. 



Table 1. Compiled catalogue of l.o.s. candidates of positive mass- 
light correlations with projected l.o.s. significance higher than Scr 
in the lensing map (shown as S/N^); S/NpM lists the correspond- 
ing signal-to-noise in the product map. Lines of sight with no 
product map match are not listed. The equatorial coordinates 
(J2000.0) arc the peak positions in the product map, redshift 
constraints arc estimates with 68% credibility from a thorough 
l.o.s. modelling in Sect. 3.2. A dagger! denotes redshift estimates 
from a two component model rather than a one-component mass 
model. 
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0.68: 
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+0.10 
-0.08 
+0.15 
-0.12 
+0.35 
-0.30 
+0.20 
-0.15 
+0.40 
-0.30 
+0.25 
-0.30 
+0.25 
-0.15 
+0.30 
-0.05 
+0.30 
-0.45 
+0.25 
-0.12 



no significant (> 3.5cr) lensing mass overdensity matched by 
a stellar mass underdensity in the entire 3D map ("dark 
clump" ) . 



3.2 Line of sight modelling 

Using the product maps in Fig. 8, we identify interesting 
l.o.s. in the map to perform a more detailed analysis of the 
3D structure that takes fully into account the radial p.s.f. of 
the Wiener map. For each l.o.s., listed in Table 1, we model 
the matter distribution by assuming that it consists of one 
individual non-zero mass halo oc 1 -f 5o on the lens planes. 
The model profile in the convolved 3D Wiener map is given 
by 



Sn-i{i, So) = Sop 



(0 



(37) 



where (5o > — 1 is the density contrast of the non-zero pixel 
on the ith lens plane; all other pixels have 5^' ~ which 
amounts to average density in the corresponding epoche. 
The radial matter density dm in direction G of the Wiener 
map with noise covariance N enables us to write down the 
posterior likelihood of the mass pixel lens plane index i and 
its contrast Sq: 



P(j,5o|<5m) ocH(5o + l) X 



(38) 



exp 



i [dm -5n,(j,5o)]'N-i 



Sm(i, So] 



where the Heaviside function prior H{x) asserts a positive 
density 1 + So 0. 

Marginalising over So, we employ the MCMC technique 
to construct a probability distribution P(i|<5m) for each lens 
plane i, as shown as black lines in Fig. 9. This probabil- 
ity distribution gives the true probability of a density peak 
positioned along each l.o.s. at a given redshift. This result 
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redshift 



Figure 10. Probability distribution for rcdshifts in a two- 
component model in the absence of any signal (null distribution). 
The distributions for z\ and Z2 are not top-hats due to the a-priori 
condition 21 < 22. 



correctly accounts for the broad radial p.s.f that has ham- 
pered the interpretation of the maps presented up to this 
point. The blue error bars denote a 68% confidence regions 
about the median redshift of the non-zero pixel. For com- 
parison, the density contrast from the M, map for the same 
pencil beam is shown as grey filled bars; the beams have a 
radius of ~ 40 arcsec. The bars are normalised such that the 
maximum density reaches the half height of the diagram. 

Motivated by the results of the product map in Fig. 8, 
we also test an alternative two-halo mass model for most 
of the 1.0. s. , which seem to require a more complex mod- 
elling when comparing the one-component p.d.f. to the M« 
distribution in the pencil beam. The posterior likelihood of 
the two-component model is similar to Eq. 38, but now con- 
tains two positions i\ and 12 and two densities. As additional 
prior, we demand that ii < 12. The resulting marginalised 
probability distributions are shown as black solid and dashed 
lines, the two error bars denote the constraints on redshifts 
in the two-component model. The orange solid line shows 
the probability distribution for the one-component model. 
Note that in a two-component model fit with no signal in 
the data, the probability densities of both components are 
not simply fiat: Due to the adopted prior zi < Z2, the null 
signal distribution has a saw-tooth shape as in Fig. 10, which 
is produced due to the correlation of the probability distri- 
butions of 21 and 2:2- 

In Fig. 9 we see good agreement between the peaks in 
the density redshift probability distribution and the over- 
densities in stellar mass. We also see clear evidence for l.o.s. 
projections in the stellar mass distribution which have a 
significant probability in the 3D density map. In almost 
all two-component model is not clearly preferred 

over a one-component model with the slight exception of 
A901q that prefers a two-component model (d.o.f. = 18) with 
Ax'^ = 1.90 (~ 68% confidence) as compared to a one- 
component model (d.o.f.=20). The best-fit one-component 
model of SWc has = 17.23 as opposed to x^ ~ 16.70 for 
the two-component model, which is statistically insignifi- 
cant. For SWc, however, a model with all mass concentrated 
a,t z = 0.165, the distance of SWa/b, is excluded with ~ 68% 
confidence (x^ = 18.69; d.o.f.=21). Combining this infor- 
mation with our knowledge of the galaxy distribution, we 
favour the two component model which reveals a new high 



redshift structure SWc at z ^ 0.7, lying in projection with 
the foreground SW group at z — 0.165. 



4 DISCUSSION AND CONCLUSIONS 

This paper uses 3D lensing data from the STAGES project 
to study the spatial distribution of matter density and galax- 
ies in the Abell 901/2 supercluster field. We have presented 
a non-parametric lensing-based 3D mass density reconstruc- 
tion. This application of the 3D mapping algorithm to data 
is the first of its kind as earlier attempts (Massey et al. 2007; 
Taylor et al. 2004) to do cosmography with 3D lensing data 
focused on a reconstruction of the gravitational potential 
rather than the mass density itself. In addition to directly 
mapping the 3D density field, the new algorithm presented 
in this paper accounts for redshift bias, a term previously ne- 
glected in the literature, and it optimally includes all lensed 
sources including those without reliable photometric redshift 
estimates.^ 

We present both our raw Wiener reconstruction of the 
STAGES 3D density field in Fig. 7 and a cleaned visuali- 
sation of the 3D field in Fig. 6. Using the visualisation is a 
valid tool to describe the structure in the field, in particular 
because, in this instance, the z-shift bias has been correctly 
removed. Further it is this cleaned map which should be 
compared to previous 3D gravitational potential maps in 
the literature where the displayed isophotal potential lev- 
els have been optimised to highlight structures. However we 
urge caution in using any of these cleaned visualisations for 
scientific analysis as they are very noisy estimators of the 
density field and do not show the redshift error in the map or 
account for the broad radial smoothing of the field, shown in 
Fig. 1. However we have created a striking visualisation tool 
of the 3D dark matter Universe whose impact should not 
be underestimated even if its direct use for scientific anal- 
yses is somewhat limited. Note, as discussed in STH09, for 
deeper observations and higher number densities of sources, 
the radial smearing effect is reduced but it will always be 
present as a generic feature of this Wiener filter 3D lensing 
method, no matter how good a data set is analysed. A fu- 
ture improvement, albeit not explored for this paper, may 
be to regularise the Wiener filter, e.g. by a prior that forces 
structure to be confined to one lens plane. Although this 
resulted in a more realistic estimate for the 3D matter dis- 
tribution, it would not change the considerable uncertainties 
in the visualisation. 

From the raw 3D Wiener lensing reconstruction we in- 
fer that most of the structure in the STAGES field belongs 
to the Abell supercluster at z ~ 0.165 since most l.o.s. den- 
sity profiles in the map are locally peaked at around this 
redshift. There are a few exceptions however, of which the 
most prominent is in the SW group region. Here, we find a 

^ In principle, the 2-shift bias is also be present in previous 
Wiener filter maps of the 3D gravitational potential. That it was 
not recognised in the method of Taylor et al. (2004) may be due 
to the difference in the Wiener filter, which processes all lines 
of sight individually, ignoring correlations of statistical errors be- 
tween distinct lines of sight. This is suboptimal in terms of signal- 
to-noise in the resulting map, but on the other hand might result 
in a less severe 2-shift bias. 
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Figure 9. Probability distribution of redshifts in a one-component (only one black solid line or orange solid line) or a two-component 
(solid and dashed black line) mass model for the radial density profile into different directions of STAGES. Grey bars indicate the stellar 
mass Mt within a pencil beam of ~ 40arcsec radius in units of Mq/i"^ (numbers on right edge). The blue error bars denote a 68% 
confidence regions about the median redshift of the mass concentration. 
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strong peak at around 2 ~ 0.7 apparent in the raw Wiener 
map. 

In order to interpret the Wiener map result we have 
compared the total mass to the stellar mass in the field us- 
ing a product map technique outlined in Sect. 2.8. In the 
resulting product map in Fig. 8, significant matches or mis- 
matches of density fiuctuations between total mass and M* 
from the selected redshift range are identified and followed 
up with a more detailed analysis in Sect. 3.2. The conclu- 
sions we can draw from this analysis are as follows: 

(i) We find no detection of "dark clumps" in the 3D map, 
i.e. matches of lensing mass overdensities with stellar mass 
underdensities, which have been discussed in the literature 
(von der Linden et al. 2006; Erben et al. 2000; Massey et al. 
2007). 

(ii) As expected from Fig. 3, A901b is well confined in 
redshift by the 3D lensing data, and the estimated matter 
distribution in the pencil beam fits well to the M, distribu- 
tion. 

(iii) A901a, the second massive cluster in the field is less 
well constrained in redshift, possibly explained by the pres- 
ence of an extra mass peak at around 2: '-^ 0.3 as seen in M* 
and also highlighted in the product map (Fig. 8, top right 
panel). 

(iv) The SW group pencil beam is equally well fit by a 
one or two-halo model in projection. However the single-halo 
model exhibits an additional peak at higher redshift and is 
slightly inconsistent with all mass concentrated at z = 0.165. 
We therefore conclude that along this l.o.s. there are indeed 
two structures; the well documented foreground SW group 
at z — 0.165 and a second possible mass concentration at 
z = 0.7 ± 0.3, which is confirmed by a corresponding stellar 
mass peak. SWc is also the most prominent high- 2; peak in 
the raw Wiener map. This finding was missed in the past 2D 
lensing analyses of Gray et al. (2002); Heymans et al. (2008) 
and the smoother 3D map of the gravitational potential in 
Taylor et al. (2004). 

(v) We find a similar result for the the A901a X-ray group 
known to be infalling on A901a at a redshift of 2 ~ 0.17. 
In this case the one-halo model again rules out a low- 
redshift structure at moderate significance, suggesting that 
the equally well fit two-halo model best represents the data 
with widely scattered matter behind 2 — 0.165 as indicated 
by the M, map. 

(vi) A902 is equally well fit by a low-redshift one-halo 
model consistent with the known low-redshift of the cluster, 
or a two-halo projection model with one at 2 ~ 0.2 and one 
at 2 ~ 0.8. 

(vii) CBI, clearly visible as M* peak at 2 = 0.46, is weakly 
detected in the mass map, which indicates a mass feature at 
a consistent redshift 2 = 0.68lQ'4g. 

(viii) We also find evidence for mass concentrations be- 
hind the A901b cluster termed as A901b/2 and A901b/3. 
The latter is more prominent in the product map, concen- 
trated at 2 ~ 0.7. 

Note that none of these conclusions could be drawn from 
the purely 2D analysis of Heymans et al. (2008). It is the 
non-parametric, purely lensing-based modelling of selected 
l.o.s.. Fig. 9, that has to be used to investigate these signif- 
icant features. In contrast to our 3D Wiener map and com- 
parable visualisations (Hu & Keeton 2002; Bacon & Tay- 



Table 2. Constraints on the NFW virial masses M200 and the 
virial radii r2oo for A901b and SWa,b in a simple model fit 
to the 3D lensing data and additionally A901b/2-3 and SWc 
in an updated model with more components. The redshift of 
the foreground mass peaks is fixed with z = 0.165 and with 
z = 0.8,0.75,0.75 for A901b/2-3 and SWc, respectively. Within 
the fit an uncertainty of Az = 0.05 is factored in. Quoted are 68% 
credibility regions about the mean of the posterior likelihood. As 
fiducial cosmology here Qm = 0.3 and Ha = 0.7 were assumed. 

Simple model Updated model 

Peak r200 M200 »'200 M200 

[h-ikpc] [h-ilO^Mg] [h-ikpc] [^-IIOI^Mq] 



A901b 

A901b/2 

A901b/3 

SWa 
SWb 
SWc 



1140"* 



16.9 



-H5.4 
4.6 



507 
496 



+211 

228 
-1-204 
215 



2.2 
2.1 



+3.1 
1.6 
+3.0 



Qfjr + 182 

48ll^^t 
486lJ?« 



+5.0 
4.4 
+4.8 
2.3 

^■^-1.8 



14.6 
3.3 



2.0 
1.9 
1.8 



+2.9 
1.4 
+2.7 
1.3 
3.8 
1.4 



lor 2003; Taylor et al. 2004; Massey et al. 2007; Vander- 
Plas et al. 2011), information on the redshift uncertainty of 
lensing mass peaks are directly visible in this Bayesian ap- 
proach. From here we can see that even for a survey as good 
as STAGES radial distances of objects at higher redshifts 
are usually weakly constrained, most strikingly visible for 
the background objects in a two-component fit. This con- 
firms the theoretical estimate in Fig. 3, which is based on 
SIS matter haloes. Nevertheless, the influence of background 
matter on the A901b, A902 and especially SW group lensing 
signal is visible in Fig. 9; compare the foreground p.d.f. in 
the two-component model to the p.d.f. of the one-component 
model. This influence is hidden in a 2D lensing analysis. 

All in all the modelling and its comparison to the stel- 
lar mass distribution confirms that the cleaned visualisation 
of the 3D lensing map correctly identifies previously known 
structure (e.g. Heymans et al. 2008) at low redshift, but also 
seems to identify hitherto unknown structure at larger red- 
shift, most notably A901b/3 and SWc. These overlapping 
mass concentrations are a possible source of error when it 
comes to mass determinations by lensing, a method that nat- 
urally measures the integrated mass within cylinders along 
the l.o.s. To gauge this effect, we fit Navarrow-Frenk- White 
(NFW) matter density halo profiles (Navarro et al. 1997; 
Wright & Brainerd 2000) with fixed redshifts and positions 
on the sky to the 3D lensing data. A method similar to 
the one in Heymans et al. (2008) is employed and hence 
not detailed here; the redshift p.d.f. of our shear catalogue 
subsamples are utilised, and constrains from all subsamples 
are combined. In one fit assuming a simpler mass model, 
we include SWa/b, A901b and A901a (all at 2 = 0.165) 
only, and in a second fit, invoking a more complex mass 
model based on the previous conclusions, we also include 
SWc (2 = 0.75) and A901b/2-3 (2 = 0.8,0.75) as possi- 
ble carrier of mass. The fit results can be found in Table 
2. Comparing the model fits with and without 3D structure 
in projection shows a slight drop in the mass estimates for 
A901b and SWa/b but this decrease is not statistically dis- 
tinguishable from the simple mass model. This demonstrates 
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that previous 2D mass estimates in the hterature may have 
overestimated mass measures but not significantly within 
the statistical errors. This result is expected based on Fig. 9 
where single components are as equally well fit as a multiple 
component model along these lines of sight. 

There is also the possibility that the observed features 
in the background are actually belonging to the foreground 
supercluster and are moved to higher redshift in the Wiener 
map owing to the noise in the lensing signal; their mass is 
relatively small (Table 2) in the light of Fig. 3. In this case, 
the match with stellar mass at higher redshift would be a 
coincidence. 

We have presented a method that can be used to ac- 
count for the radial smearing of the 3D mass map and pro- 
duce full redshift probability distributions along individual 
l.o.s.. In future work an alternative, potentially more fruit- 
ful method could be to first identify significant lensing de- 
tections in a 2D map and then produce smaller scale 3D 
reconstructions in these interesting regions. The benefit of 
this method would be that each 3D reconstruction could be 
optimised for each source by tuning the a parameter in the 
Wiener reconstruction (see Sect. 2.2) depending of the mass 
and redshift of the source. In addition, information from 
measurements of higher order weak lensing flexion distor- 
tions could be included in the analysis to boost the signal to 
noise in the reconstruction (Velander et al. 2011). Moreover, 
for deeper surveys, which have reliable photometric redshifts 
down to 2 ~ 2 or more, we also anticipate an moderately 
improved mass density recovery at higher redshift based on 
simulations (Fig. 13 in STHOgj. 

Throughout this paper, we have neglected errors caused 
by the intrinsic alignment of galaxies, which for this analysis 
would be dominated by the 'GI' term, the anti-correlation 
of sheared background sources with the intrinsic shape of 
foreground galaxies. Based on analysis of simulations and 
data, the systematic error contributed by this astrophysical 
effect is expected to be low on average (Joachimi & Bridle 
2010) and probably negligible for the results presented in 
this paper. That said, the intrinsic alignment of galaxies in 
cluster environments has yet to be studied in data but has 
been explored using analytical prescriptions (Schneider & 
Bridle 2010). If cluster galaxies strongly align as they infall 
this could result in a strong GI term that could effect our 
results by lowering the lensing signal measured behind the 
foreground cluster. The opposite however could also occur 
whereby galaxy merging events in the dense cluster environ- 
ment randomise the intrinsic ellipticities resulting in a very 
weak GI term. As discussed in STII09, the contribution of 
GI and II correlations to the shear field can in principle be 
incorporated into the minimum variance estimator for the 
3D map. At this stage, however, we expect this effect to be 
small. This uncertainty is something that will be investi- 
gated further using the STAGES data in future work. 

The driving motivation behind the development of our 
3D reconstruction technique was to enable an unbiased 3D 
comparison of mass and light. Dark haloes for example 
would only be detected in this manner. However with the 
detailed analysis of noise and the radial p.s.f in the 3D lens- 
ing reconstructions presented in this paper, we have demon- 
strated how inherently noisy the process is. Even with one 
of the best HST data sets and a very massive foreground 
system we have found it difficult to reconstruct the 3D mass 



distribution at any high degree of spatial and redshift reso- 
lution using weak lensing alone, especially at higher redshift. 
Given the limitations of the method to resolve only the most 
massive structures in 3D the future direction for the appli- 
cation of this method should be to reconstruct larger scale 
structures in the 3D density field using more heavily spa- 
tially smoothed data. With upcoming wide area surveys such 
as CFHTLenS, DES and VST-KIDS, which are several or- 
ders of magnitude wider than the STAGES survey presented 
in this paper, we can expect higher quality 3D resolution re- 
constructions on degree scales because on these scales the 
significance of modes in a 3D mass density reconstruction 
are increased (STH09). 



ACKNOWLEDGEMENTS 

We thank Peter Schneider and Jan Hartlap for useful dis- 
cussions. We also thank the extreme snowy weather condi- 
tions that kept PS and CH snowed in at the Royal Observa- 
tory, forcing them to finally finish this paper. PS acknowl- 
edges supported by the European DUEL Research- Training 
Network (MRTN-CT-2006-036133) and by the Deutsche 
Forschungsgemeinschaft under the project SCHN 342/7- 
1. CH acknowledges support from the European Research 
Council under the European Communities Seventh Frame- 
work Programme, ERG grant agreement number 204185. 
TS acknowledges support from the Netherlands Organisa- 
tion for Scientific Research (NWO). MG was supported by 
an STFC Advanced Fellowship and also acknowledges fund- 
ing by the STAGES NASA grant. KJ is supported by the 
Emmy Noether-Programme of the German science founda- 
tion DFG. BH is grateful for support from the Science and 
Technology Facilities Council (STFC). Support for STAGES 
was provided by NASA through GO- 10395 from STScI op- 
erated by AURA under NAS5-26555. 



APPENDIX A: THE STAGES LENSING 
ANALYSIS 

This Appendix documents the steps taken to confirm the 
reliability of the weak lensing shear catalogue that forms 
the basis of the 3D reconstruction presented in Sect. 3. 

Al Shear catalogue: Sample I and Sample II 

The HST STAGES data were analysed independently by 
two authors (CH and TS). This involved fully independent 
data reduction of the Advanced Camera for Surveys (ACS) 
observations, independent source extraction, p.s.f. modelling 
and independent shear measurement pipelines based on the 
Kaiser et al. (1995) (KSB) method. The full HST lensing 
analysis is detailed in Heymans et al. (2008) (for CH) and 
Schrabback et al. (2007, 2010) (for TS) with both KSB- 
based shear measurement pipelines compared in Table Al 
of Heymans et al. (2006). 

The main differences between the two lensing analyses 
are the p.s.f. modelling method and the galaxy selection. CH 
assumes the p.s.f. temporal variation is slow over time and 
splits the HST data in 7 groups to achieve good temporal 
sampling of the p.s.f.. The resulting semi-time dependent 
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p.s.f. model is then applied to the data. CH also chooses to 
maximise the signal-to-noise in the lensing catalogue fitting 
the isotropic correction term P'' as a function of galaxy size 
(see Heymans et al. (2008) for more details). 

In the analysis of the dense stellar fields, TS has shown 
the HST ACS p.s.f. can vary rapidly with time. To account 
for this, TS employs a fully time-dependent correction for 
the spatially and temporally varying p.s.f. based on indi- 
vidual exposures, using a p.s.f. template library obtained 
from stellar fields (Schrabback et al. 2007). In addition, 
TS applies corrections for charge transfer inefficiency and 
signal-to-noise dependent shear calibration bias as detailed 
in Schrabback et al. (2010). In this analysis of the STAGES 
data TS employs conservative galaxy selection criteria in- 
cluding a high cut on the shape measurement signal-to-noise 
ratio (Erben et al. 2001) S/N > 2.7. 

The resulting catalogues differ in the number density 
of objects and the noise on the shear measurement quan- 
tified through the root-mean square variation of the mea- 
sured galaxy ellipticities. CH, hereafter Sample I, contains 
65 galaxies per square arcmin with = 0.26. TS, here- 
after Sample II contains 34 galaxies per square arcmin with 
a, = 0.29. 

A2 Photometric redshifts: COMBO-17 sample 

The STAGES HST data is complemented by the 17-band 
optical COMBO-17 survey. Using a combination of broad 
and narrow-band filters, COMBO-17 uses template-fitting 
to estimate photometric redshifts (Wolf et al. 2004). For the 
A901/902 field the redshift error in units of (Jz/{'i + z) is 
estimated at « 0.02 to i? < 23 with an increase to 0.035 
in the interval R = [23,24] (see Gray et al. 2009). Outlier 
rates also increase with magnitude: we estimate roughly 1% 
outliers with |A2|/(1 + z) > 0.10 at i? < 22 and 7% outliers 
with |Az|/(l + z) > 0.15 aX R = [23, 24]. Above a magnitude 
f? = 24 we do not use COMBO-17 redshift estimates due to 
rapidly increasing errors and outlier rates. 

Owing to the width of the lensing kernel in redshift 
space, the accuracy of the COMBO-17 redshifts are more 
than sufficient. Redshift estimates, however, only exist for a 
small fraction of the galaxies (13 galaxies per square arcmin 
for Sample I and 10 galaxies per square arcmin for Sam- 
ple II). For the remaining galaxies we assign a magnitude 
dependent redshift distribution given by Schrabback et al. 
(2007) (see Fig. Al). All this information is combined for 
the final mass map. 

The stellar mass estimates used in this analysis derive 
from low resolution 17-band spectra fits to parameterised 
star formation history models as described in Borch et al. 
(2006). 

A3 Quality control; comparison of E and B mode 
signals to pure noise maps 

A reconstruction of signal and noise, based on the shear 
Sample I and shear Sample II, is carried out following STH09 
with the resulting map is shown in Fig. 6 for Sample I and 
Fig. A2 for Sample II. 

We find that both samples detect the same significant 
structures at the positions of the known galaxy clusters 
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Figure Al. Distribution of sources in subsamples as used in the 
3D reconstruction. The solid line shows the distribution of pho- 
tometric redshifts for sources with photo-z's, the other lines are 
estimated distributions for three magnitude bins (Sample I; Sam- 
ple II is only marginally different) including all other sources: 
23 ^ mp606W < 25 (dashed), 25 ^ mpeoew < 26 (dashed- 
dotted) and 26 ^5 "ipeoew < 27.5 (dotted). Sources with pho- 
tometric redshifts are further sub-divided into 21 photo-z bins. 
Distributions have arbitrary normalisations. 

A901a, A901b and A902, in addition to a detection of a 
group of galaxies called the South West group (SW) as seen 
before in earlier 2D reconstructions of the projected mat- 
ter density (Heymans et al. 2008; Gray et al. 2002). As ex- 
pected from the galaxy number counts in the two samples, 
the signal-to-noise in the maps is higher for Sample I than 
Sample II. All structures that differ between the maps are 
at a level consistent with noise. 

For an analysis of the systematics, we count the num- 
ber of peaks and troughs in the E- and B-mode S/N 3D 
maps given signal-to-noise range and compare this to the 
expected numbers in pure noise maps. The noise statistics 
are obtained from the randomised data. In order to have an 
easily computable and objective statistic for "peaks" and 
"troughs" , we count the number of grid pixels in the Wiener 
map under the condition that the pixels are local maocima 
within a cuboid volume centred on the pixel. The size of 
the cuboid is ±2 pixels in transverse and ±2 lens-planes in 
radial direction. For S/N-thresholds > 3o- this definition fits 
reasonably well the number of peaks one would select by eye. 
Troughs are considered negative peaks (local minima with 
negative signal-to-noise) . 

Fig. A3 shows the results for Sample I and Sample II. 
Clearly, for thresholds > 3.5(t, E-mode peaks (troughs) are 
above the average noise levels and above Itr-noise variance 
for > 4.0(7, whereas B-mode counts are consistent (Sample 
I) with the noise expectation. From the noise realisations, we 
estimate the covariance of the noise counts which is used for 
a x^-test, peak and trough counts combined, to check con- 
sistency with noise. For Sample I we find x^/d.o.f. = 177.4 
(E-mode map) and x^/d.o.f. — 0.40 (B-mode map), show- 
ing that the E-mode map is inconsistent with noise, while 
the B-mode map is consistent with noise. Corresponding 
values for Sample II are /A.o.i. = 220.0 (E-mode) and 
X^/d.o.f. = 0.29 (B-mode) leading us to a similar conclu- 
sion. 

We therefore find the B-modes for both Samples to be 
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Figure A2. Comparison of the mass reconstruction from SAM- 
PLE I and SAMPLE IL On-sky projection of Wiener map using 
SAMPLE II (intensities in the background are maximum density 
contrast along l.o.s.). Black solid lines are signal-to-noise levels 
2.5(7, 3(7, 4(7, 5(7 and 6(7 of the E-mode map. The blue contours 
show the signal-to-noise in the B-mode map. For comparison, the 
magenta contours encircle B-mode signal-to-noise peaks in the 
SAMPLE I map. The values on the a;-axis denote the right as- 
cension in degrees and the declination for the y-axis (J2000.0). 
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Figure A3. Number of peaks and troughs within a given signal- 
to-noise bin {top panel: Sample I, bottom panel: Sample II). The 
orange bar in the middle shows the mean number of noise peaks 
and their r.m.s. variance as determined from noise realisations 
(10^), solid bars are statistics of the E-mode map, whereas shad- 
owed bars show the statistics of the B-mode map. The exact def- 
initions of pixel peaks and troughs is given in the text. Note that 
the noise statistics of peaks and troughs is within numerical errors 
identical so that both are represented by the same orange bar. 



low and consistent with noise. This is potentially a surprise 
for Sample I, bearing in mind the rapid variations in p.s.f. 
ellipticity seen in Schrabback et al. (2007) but unaccounted 
for in the Sample I analysis. The success of the semi-time 
dependent p.s.f. subtraction method of Sample I may be ex- 
plained by the fact that HST was significantly under-focus 
during the STAGES observations with the result of a rela- 
tively stable p.s.f., once averaged over intra-orbit variations. 
Owing to the larger Sample I, there are peaks up to 7(t, 
unseen for Sample II. 

Note that in the high-S/N regime the number of peaks 
or troughs becomes rare and statistically Poissonian (one 
or no peak) which may not be properly reflected by the 
symmetric error bars in the figure. 

A4 Redshift Scaling 

Recent lensing analysis of ground-based data using the KSB 
method found errors in the expected redshift scaling of the 
measured cosmic shear signal (Kilbinger et al. 2009). The 
reason for this bias is unknown but is thought to lie in a weak 
signal-to-noise and galaxy-to-psf size dependent calibration 
bias inherent to shear measurement methods (Massey et al. 
2007; Bridle et al. 2010; Semboloni et al. 2009). 

As a further quality check we test the redshift scaling 
of the KSB-measured space-based lensing signal for Sample 
I and II. We select a set of peaks that are known to be at 
redshift z ^ 0.165 (A901a/b/Q?, A902/1 and SW/a/b). For 
each peak the mean tangential shear within an aperture, of 
radius between one and two arcmin, is measured from the 
data: 



(Al) 



where 0^ is the l.o.s. direction of the aperture centre, Oi 
the angular position of a source in the common complex 
notation and ei the complex ellipticity of the source. Only 
sources within the aperture and within some redshift range 
are considered. The statistical uncertainty of the estimator 
for the mean shear is the variance of the mean as obtained 
from the sources included in the sum. The mean redshift, z, 
of included sources determines the redshift associated with 
7t. The statistics from all lensing peaks is combined and 
weighed by their statistical uncertainties. Thus we are stack- 
ing the cluster signals to increase the signal-to-noise of this 
test. 

Fig. A4 shows aperture statistics for Sample I and Sam- 
ple II by plotting the real parts of 7t (E-mode) and the imag- 
inary parts (B-mode; radial shear) separately. If the lensing 
signal within the apertures stems in essence only from struc- 
ture located at the cluster redshift z ~ 0.165, we expect a 
redshift scaling dictated by the lensing efficiency (with fixed 
lens distance): 



lt{z) (X 



/k(Ws 



/k(lfs) 



/k((Ws) - Wh) 
/k(K)) 



(A2) 



where f]^{w) is the angular diameter distance as function of 
radial comoving distance w; w\-i, here constant, and Ws are 
the halo and average source distance, respectively. The last 
step is only approximately correct. This approximation is 
used for the fit in Fig. A4. 

This behaviour is independent of the intrinsic density 
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Figure A4. Test of redshift-scaling of the lensing signal in the 
data. Mean tangential shear as function of source redshift within 
an aperture of radius 2 arcmin (A901a and A901b) or 1.4 arcmin 
(A901a, A902/1 and SWa/b) centred on the peak centres; the 
signal for all four apertures is averaged. The sources are binned 
in redshift for the samples with photometric galaxy redshifts (z ^ 
1.2). The three highest 2-bins belong to the magnitude samples 
for which redshifts correspond to the mean as determined from 
the redshift distribution of sources. With most matter inside the 
aperture located at the same redshift, the mean tangential shear 
(solid points) should scale according to the black solid line in 
the plot. The shaded region around the line is the Icr-confidence 
region of the amplitude. Open crosses denote the average radial 
shear inside the aperture, which should be vanishing. The shaded 
red area marks the Icr-area for the mean B-mode over all redshifts. 
Top panel: Sample I, bottom panel: Sample II. 



profiles of the clusters, albeit its ampfitude depends on the 
cluster details. As we are solely interested in verifying the 
correct redshift scaling behaviour, we fit Eq. A2 with the 
amplitude as unknown variable to the E-mode data points 
which for xVd-o.f. = 0.46(0.50) is an excellent fit to Sample 
I (Sample II). 

B-modes should be consistent with zero for all redshifts. 
This we test by fitting a straight horizontal line, with un- 
known amplitude, to the average radial shear. We find that 
the fit is both consistent with zero on a Icr-level and well 
fitted by a model with redshift-independent signal, yielding 
xVd-o.f. = 0.97(0.87). 



A5 Conclusion 

In conclusion, we find both shear catalogues to have a low 
level of systematics in the B-mode channel consistent with 
noise and a redshift-scaling of the shear signal that is consis- 
tent with ACDM cosmology. Moreover, we find no significant 
differences between the E-mode channel maps of Sample I 
and Sample II. As Sample I yields an overall more signifi- 
cant map owing to the higher number density of sources, we 
present results in the main section of the paper only from 
the data of Sample I (CH). 
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